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Abstract 



We consider how well current theories can predict neutrino mass and mixing param- 
eters, and construct a statistical discriminator which allows us to compare different 
models to each other. As an example we consider see-saw models based on fam- 
■ ily symmetry, and single right-handed neutrino dominance, and compare them to 

I> I each other and to the case of neutrino anarchy with random entries in the neu- 

trino Yukawa and Majorana mass matrices. The predictions depend crucially on 
the range of the undetermined coefficients over which we scan, and we speculate 
on how future theories might lead to more precise predictions for the coefficients 
' and hence for neutrino observables. Our results indicate how accurately neutrino 

masses and mixing angles need to be measured by future experiments in order to 
discriminate between current models. 



^ . 1 Neutrino oscillation physics is one of the most rapidly developing areas of 

particle physics [1]. Current atmospheric neutrino data are consistent with 
I t'/x ^ i^T oscillations [2] with sin^ 26'23 > 0.88 and 1.5 x 10~^ eV^ < |Am32| < 

^ ; 5 X 10"=^ eV2 (90% CL) [3], where Am^. = - in units of eV^. CHOOZ 

data [4] limits sin^ 26113 < 0.1 - 0.3 over the SuperKamiokande (SK) pre- 
ferred range of Am32. The solar neutrino problem can be solved by four dif- 
ferent combinations of oscillation parameters, three of which are based on the 
Mikheyev-Smirnov-Wolfenstein (MSW) mechanism [5], while the fourth in- 
volves neutrino oscillations in vacuum (VAC). These are the well-known large 
mixing angle (LMA), the small mixing angle (SMA) solutions and a third so- 
lution involving a large angle but a lower Am^2 (LOW solution). The best-fit 
points of the various solutions are approximately {Aml2, sin^ 2^i2,name) [3]: 
(3.3 X 10~^0.78, LMA MSW); (5.1 x 10"^, 0.0027, SMA MSW); (10-^0.93, 
LOW MSW); (8 x 10"^°, 0.93, VAC). 

Recently the SNO collaboration [6] has presented its first data on the charged 
current flux of the solar neutrinos, which supports both the hypothesis 
of active neutrino oscillations and the standard solar model [7]. Furthermore 
results of recent global analyses, performed just days after the SNO result 
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[8] show that the most favoured solutions are either the LMA or the LOW 
solutions, with a slight preference for the LMA solution. The SMA solution is 
currently disfavoured at about 3 cr, while the VAC regions are disfavoured at 
various confidence levels. 

Because oscillation experiments do not measure absolute values of neutrino 
masses, in the numerical part of this work we will investigate the following 
quantities, 



R = |Am2i|/|Am^ 
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If U is regarded as the lepton mixing matrix, then, in the limit ^ 1, to 
very good approximation sc ~ sin^ 26*13, =^atm ~ sin^ 26*23, -s© ~ sin^2^i2. In 
the following we will call ^13, 6^23, ^12 the "CHOOZ", atmospheric, solar angle, 
although strictly speaking the various experiments measure sc, Satm, Sq in 
the hmit sc < 1- Note that for the LMA MSW solution R is of the order of 
0(10-2). 

Following the experimental confirmation of the atmospheric neutrino problem 
by the SK collaboration [2], there has been a fiood of papers on models of 
neutrino masses [9]. In order to plan future neutrino experiments it would be 
very helpful to know how accurate future measurements would need to be to 
usefully discriminate between the different models, and eventually to decide 
which (if any) of all these theoretical attempts is the correct one. In order to 
address this question it is important to understand how well current theories 
predict neutrino masses and mixing angles. In view of this we find it rather 
surprising that although hundreds of different models have been proposed 
there has so far been no attempt to evaluate how well the different models 
actually predict the neutrino parameters. The main purpose of this paper is 
to address this issue. 

In most of the neutrino mass models [9] based on the seesaw mechanism [10] 
one usually assumes that there is some underlying broken symmetry which 
generates a small dimensionless expansion parameter A. The entries in the neu- 
trino mass matrix then will depend on A to some - presumably - high power. 
Such models are called "texture" models, since they produce approximate ze- 
roes in various entries in the mass matrices under consideration (textures). 
Coefficients in front of the expansion parameter, which are not predicted in 
these models, are then usually assumed to be order 0{1) numbers, and their 
infiuence on the final result is neglected. In this work we are going to inves- 
tigate the relative importance of the unknown order C(l) coefficients for the 
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first time systematically. 



Our basic idea is to construct the following simple neutrino mass model dis- 
criminator which we define qualitatively in the following sense: calculate the 
leading structure (given as powers of A which are fixed by symmetry) of any 
interesting model and scan "randomly" over the unknown parameters. From 
these samples one can construct distributions of measurable neutrino parame- 
ters, which can be compared to experiment. A model is then considered "good" 
if the peaks of the distributions coincide with (or are close to) the experimen- 
tally preferred value. Such a discriminator is by its nature qualitative rather 
then quantitative, and it is unclear how much weight should be attributed to a 
model that predicts all quantities close to the experimentally preferred values. 
The main problem of the approach is that we have treated all the coefficients 
on an equal footing, and one can imagine some extra condition which could 
become known in the future which could single out a particular value for a 
"random" parameter, which would lead to predictions in the tails of the dis- 
tributions, rather than near the peaks. However, we would like to argue that 
although a priori such a possibility cannot be strictly excluded, it might be 
considered to be contrary to the spirit of texture models. 

In our analysis we only consider "high energy scale" models based on the 
seesaw mechanism [10], which are otherwise hard to test. From the "testabil- 
ity" point of view, low-energy scale models of neutrino mass (see, for example 
[11,12]) might be considered preferable. However, even if one of the low-energy 
approaches is the true neutrino mass model, a confirmation certainly lies a 
number of years in the future, and so our methods of discrimination could be 
useful even in that case. 

2 For definiteness we consider [13] models based on a C/(l) family symmetry. 

The idea of such a symmetry is that the three families of leptons are assigned 
different U{1) charges, and these different charges then control the degree of 
suppression of the operators responsible for the Yukawa couplings, leading to 
Yukawa matrices with a hiearchy of entries, and approximate "texture" zeroes. 
As usual it is assumed that the U{1) is broken by the VEVs of some fields 
9, 6 which are singlets under the standard model gauge group, but which have 
vector-like charges ±1 under the f/(l) family symmetry. The U{1) breaking 
scale is set by < 6* >=< 6 >. Additional exotic vector matter with mass 
My allows an expansion parameter A to be generated by a Froggatt-Nielsen 
mechanism, 

^ = ^-^-0-22 (5) 

Mv Mv 



where the numerical value of A is motivated by the size of the Cabibbo angle. 
Small Yukawa couplings are generated effectively from higher dimension non- 
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Table 1 

Flavour charges (FC) for four models, as discussed in the text. 



renormalisable operators corresponding to insertions of Q and Q fields and 
hence to powers of the expansion parameter in Eq.5. The number of powers 
of the expansion parameter is controlled by the charge of the particular 
operator. The lepton doublets, neutrino singlets, Higgs doublet and Higgs 
singlet relevant to the construction of neutrino mass matrices are assigned 
charges /j, np, /i„ = and a. From Eq.5, the heavy Majorana matrix 
elements are M^^^ ~ Aij}^''^'-^'"'^^"^ The neutrino Dirac matrix elements are 
ml^ ~ aijAl^*+"-'', where A^j and undetermined coefficients. For details 

see [14]. 

The theory does not determine the coefficients and which we take to 
be free parameters which are chosen randomly over some range. The physical 
light effective Majorana neutrino mass matrix which results from the see-saw 
mechanism [10] is given by m^L = ^LRMj^j^mj^ji, and itill is diagonalised 
by a matrix U and the eigenvalues are the light neutrino masses. In the 
basis in which the charged lepton mass matrix is diagonal, U is just the lepton 
mixing matrix. 

In Table 1 we give four examples of models based on different choices of 
flavour charges. Models FC1-FC3 arc examples of single right handed neutrino 
dominance (SRHND) [14] which naturally lead to small values of R, typically 
i? ~ (see below). The essential idea of SRHND is that one of the right- 
handed neutrinos contributes dominantly to the 23 block of rriLL- This results 
in the approximate vanishing of the the 23 sub determinant of rriLL, and hence 
one expects a hierarchical spectrum of masses, m^^ <^ my^. Without SRHND 
a large mixing angle O^z ~ 7r/4 would typically result in two eigenvalues of 
the same order of magnitude m2 ~ ma, and the presence of a mass hierarchy 
could only be achieved at the expense of some fine-tuning. The amount of 
fine-tuning required depends on the value of the mass hierarchy. 

FCl and FC3 are taken from [15] where it is shown that they lead to the 
expectation of a large solar angle. FCl leads to dominant (from the third right- 
handed neutrino) and subdominant contributions to rriLL [15] (order 0{1) 
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coefficients suppressed for simpficity): 



^ A2 A^^ 
A2 1 1 
A2 1 1 



+ 



^A^ A2 A2^ 
A2 A2 A2 
A2 A^ A2 



(6) 



We constructed FC2 as an example of a flavour model predicting a small solar 
angle. FC2 gives dominant and subdominant contributions: 



ml? - 



( A^ A^ A^ 
A^ A^ A^ 
yA^ A^ A3 



+ 



^A^ A^ A^^ 
A^ A^ A^ 
A^ A^ A^ 



(7) 



FC3 gives a large solar angle, and also a larger CHOOZ angle: 
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(8) 



The model FC4 is deliberately chosen not to have SRHND, but is designed 
to give the same leading order structure as mf^^, but with subleading con- 
tributions to the 23 block which are also of the same order, leading to an 
unsupprcsscd vahie of R. We will not consider renormalization group effects 
here, since it has been shown that they only lead to corrections of several 
percent [16] which are negligible compared to the effects of the unknown co- 
efficients. 

Approximate expectations for the experimentally accessible quantities (^23, 
6'i3, O12 and R = |Am^2l/l^"^23l) be estimated from the form of the 

matrices in Eqs.6-8 as discussed in [15] and are given in Table 2. 

Finally, as a fifth example we will consider neutrino anarchy [17]. In the lan- 
guage of flavour charges anarchy simply corresponds to choosing h = I2 = h 
and rii = n2 = n^, such that the effective neutrino mass matrix shows no 
obvious structure. It is therefore completely determined by the random coef- 
ficients. 

3 As in [17] we investigate the quantities R, Sc, Satm and Sq, as defined in 
Eqs.1-4. As mentioned in the introduction these are appropriate quantities 
when considering data from oscillation experiments. Recall, in the limit sq <S 
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Table 2 

Approximate expectations for 023, 6*13, Ou and i? for the four different choices of 
flavour charges defined in the text. 

1, to very good approximation sq ~ sin^ 2^13, Satm ~ sin^ 26*23, -5© ~ sin^ 2^12. 
As a test of our numerical procedure we ran a sample calculation with one 
million randomly chosen points for all five models imposing the same cuts as 
used in [17] {R < 1/10, sc < 0.15, Satm > 0.5 and Sq > 0.5). For this test 
wc restricted the range of the coefficients to be in the interval [-1,1] in order 
to be able to compare our calculation to the results presented in [17]. Note, 
however, that this range is not the appropriate one for texture models, as 
discussed below, and was used only in this test calculation. We have checked 
that our results for the anarchy matrix are statistically consistent with the 
original anarchy calculation in [17], namely that approximately 3% of the 
points pass all the cuts. By comparison, with the same procedure applied to 
model FCl in Table 1 25% of all points pass the cuts. 

In texture models, however, it is necessary to choose the range of the coeffi- 
cients Aij and differently from [-1,1]. As explained in more detail in section 
4, where the importance of the coefficients is discussed, in order to not destroy 
the texture structure imposed by the expansion parameter, we have chosen the 
range of the coefficients to be ±[\/2A, l/\/2\]. In other words the magnitude 
of each of the coefficients is selected randomly in the range 0.663-1.508 and 
its sign is also selected randomly. With the structure preserving range only 
1.8% of the neutrino anarchy points pass the cuts, while for the case with 
structure FCl 56% of the points now pass the cuts. Structure in the neutrino 
mass matrix is clearly preferred over anarchy, since more points pass the cut 
in that case. 

In Figure 1 we give distributions for Satm, Sq,sc and R for 10*^ points and 
five different models, using random coefficients in the range ±[\/2A, 1/V2\]. 
Note, that the plots are on a log-linear scale. Using a semilogarithmic scale is 
suitable for quantities which a priori can be distributed in intervals ranging 
over several orders of magnitude, as can be expected for the observables we 
are interested in. Choosing a linear scale a model predicting Sq = 0.1, for 
example, could not be distinguished from a model predicting Sq = 10^'^. Thus 
our choice of scale and our definition of a "good" model corresponds to a rather 
mild requirement. Instead of asking for an acceptance of a given model at some 
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Fig. 1. Theoretical distributions for the predictions of neutrino mass and mixing 
parameters for five selected see-saw models: FCl (full), FC2 (dot-dashes), FC3 (thick 
dots), FC4 (thin dots), anarchy (dashes). Matrix coefficients are randomly chosen 
in the interval it [\/2A, 1 / ^/2X\ . The vertical axis in each panel (deliberately not 
labelled) represents the logarithmically binned distributions with correct relative 
normalisation for each model, with heights plotted on a linear scale in arbitrary 
units. 

(high) level of confidence, we ask for the minimum requirement that models 
should give at least the correct order-of-magnitude estimate for observables. 
Any model failing even this mild test will fail even more badly once subjected 
to a more sophisticated numerical analysis. 

The dashed line and the full line correspond to the case of neutrino anarchy 
and the FCl model as in Table 1. If we focus only on these to begin with we 
see that they both give a peak at Satm — 1, and Sq — 1, but the FCl model 
gives also a peak at Sq — 0.4, and disfavours Sq < 0.1. 

The main difference between anarchy and FCl is in the CHOOZ angle where 
anarchy again gives a peak at = 1, while the FCl model has a peak near 
the small value sc — 0.02. The R value for the FCl model also prefers smaller 
values than the case of anarchy, albeit with a broader peak, and may be 
consistent with either the LMA or the LOW solutions. Clearly more accurate 
experimental determinations of physical parameters will eventually provide 
a powerful discriminator between the case of neutrino anarchy and neutrino 
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structure. 



The remaining lines in Figure 1 correspond to the different flavour models 
FC2, FC3, FC4. All the models predict large values of Satm and so cannot be 
distinguished this way. The FC2 model (dot-dash) gives a peak at Sq — 0.01, 
and therefore prefers the SMA MSW solution, consistent with the analytic 
expectations discussed previously [15]. It also has peaks near sc = 0.01 and 
R = 10~^ consistent with the analytic expectations. The FC3 model (thick 
dots) is consistent with the LMA or LOW solutions, similar to the FCl model 
(solid), since it peaks at large Sq and the R peak is in the same region, 
albeit with a different detailed shape. However, unlike FCl, the FC3 model 
has a peak at the larger value sc = 0.3 and so may be distinguished by a 
measurement of ^13. Finally the FC4 model (thin dots) is designed to give 
the same leading order contributions to itill as FC3 (thick dots), differing 
only in the subleading contributions, so that FC4 does not possess the feature 
of SRHND. This subtle difference between the two models gives rise to a 
quantitative difference since the peak in R for the FC4 model is shifted to the 
right compared to that in the FC3 model. 

4 We emphasise that in all flavour models of the kind we consider [13] the 
coefficients Aij and not specified. In order for these models to make 
sense it is implicitly assumed that the coefficients are of order unity to within 
an error determined by the expansion parameter. The range of coefficients ap- 
plied previously, namely ±[\/2A, 1/a/2A], was chosen in order that the struc- 
ture imposed by the flavour charges and the expansion parameter A was not 
destroyed. We re-emphasise that this is not guaranteed, but a minimum re- 
quirement for texture models in order to make sense physically. We think that 
future theoretical progress in flavour models of neutrino masses lies in being 
able to predict these coefficients more accurately, and Figure 2 demonstrates 
the importance of this point. 

In Figure 2 we have plotted the distribution in the solar angle for the case 
of FC2, but with 3 different ranges of the coefficients. As the range of the 
coefficients is reduced, the central peak of the distribution becomes more pro- 
nounced, with smaller side-peaks emerging as a result of accidental cancella- 
tions due to the arbitrary signs of each element which we allow in our Monte 
Carlo. Thus Figure 2 shows the accuracy of the prediction for the solar angle 
which may be possible in some future theory where the coefficients arc speci- 
fied more accurately. Analagous results are obtained for the other observables. 

5 To conclude, we have considered how well current theories, in particular 

those based on the see-saw mechanism and SRHND, can predict neutrino 
mass and mixing parameters, and have proposed a qualitative discriminator 
based on logarithmic distributions calculated from a random scan over the 
(presently unknown) coefficients in the neutrino Yukawa matrix and the neu- 
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Fig. 2. Plots of Sq for the case of FC2 with three different ranges of the coefficients 
Aij and a,j. Full line: ±[V2X,1/V2\]. Dotted line: ±[0.82,1.18]. Dash-dotted line: 
±[0.95,1.05]. 

trino Majorana matrix, leading to the distributions in Figure 1. This paper is 
the first serious attempt to show how the uncertainty of the coefficients of the 
input see-saw matrices feeds through to the uncertainty in the predictions of 
experimental neutrino observables. 

What is the value of such distributions, and what conclusions can be drawn 
from them? Consider for example the distributions in the CHOOZ angle sq 
corresponding to the third panel of Figure 1. It is clearly seen that three 
of the models FC3 (thick dots), FC4 (thin dots) and anarchy (dashes) give 
distributions which peak at values of sc which are larger than or around the 
experimental limit sc < 0.1 — 0.3. Indeed anarchy peaks at sc — 1- This panel 
also shows qualitatively how increasingly unlikely these models will become as 
the experimental limit on Sc becomes stronger, and again anarchy fares the 
worst having the lowest number of points in the tail of the distribution. We 
should emphasise that such a plot of the distribution in sc was not shown for 
anarchy in [17]. On the other hand the models FCl (full) and FC2 (dot-dash) 
give distributions which peak at much smaller values of sc, in agreement with 
the analytic estimates in Table 2. The value of the distributions is that one can 
see pictorially how much a given experimental limit on sc favours or disfavours 
a particular model. The conclusions which one draws from such distributions 
remain a matter of taste, but at least the distributions enable such conclusions 
to be made based on the best possible information. It is manifestly clear that 
seeing a probability distribution is more informative than simply having an 
analytic estimate of the position of the peak of the distribution, and our paper 
supplies this additional information. 

Clearly from the above discussion an accurate measurement of ^13 is very im- 
portant, since it could rule out or disfavour large classes of models, and this 
may require a neutrino factory [18]. What about the other neutrino observ- 
ables? From the first panel in Figure 1 we see that an accurate determination 
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of the atmospheric angle 623 is not so important since all the considered models 
peak at maximal mixing. However, by contrast, a measurement of the solar 
angle 612 is important in resolving large classes of models, and the second 
panel in Figure 1 shows the distributions in Sq. The recent SNO indication of 
a large solar angle clearly disfavours the model FC2 (dot-dash) which peaks 
at a small solar angle according to the analytic estimate in Table 2. One may 
ask, what do we learn from the distributions that we do not already know 
from the analytic estimates of the position of the peak? The answer is that 
the distributions give us a rough idea of the uncertainties of the analytic esti- 
mates. 

The distributions in R, corresponding to the fourth panel in Figure 1, all 
have rather broad distributions which could not have been predicted by the 
analytic results, which again only predict the position of the peaks. One would 
expect that a small measured value of R would favour SRHND over other 
models where R is larger. This effect is clearly seen since the FC4 models (thin 
dots) and the anarchy model (dashes) do not have SRHND and therefore give 
distributions in R in the fourth panel in Figure 1 which are peaked towards 
larger values than the SRHND models FC1-FC3. However what could not have 
been predicted is the width of the distributions which are so broad as to make 
these distributions overlap over large regions of R. Once again the distributions 
provide this additional information which has hitherto been missing from the 
literature. Moreover, as seen from Table 2, for the models FC1-FC3 one expects 
it! to be ~ in all cases, whereas Fig. 1 demonstrates that even the position 
of the peaks in this variable can vary up to an order of magnitude. Again, the 
distributions give us an estimate of the limitations of the analytic approaches. 

We hope that the additional information provided by the distributions will 
be of use to our experimental colleagues in planning the next generation of 
neutrino experiments, but we would caution that future theories could predict 
smaller ranges of the unknown coefficients, which could demand higher exper- 
imental accuracy in the measurement of neutrino mass and mixing parameters 
than that deduced from the results presented here. Also we emphasise that we 
have treated all the coefficients as independent and uncorrelated, and some 
future theory could relate these coefficients in such a way as to give rise to pre- 
dictions in the tails of the distributions rather than near the peaks. To some 
extent this undermines the conclusions which we can draw from our current 
analysis. However we feel that the distributions which we present are a fair 
reflection of current theoretical uncertainties. 
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